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Entanglement generated between a single atom and a laser pulse 
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We quantify the entanglement generated between an atom and a laser pulse in free space. We find 
that the entanglement calculated using a simple closed-system Jaynes-Cummings Hamiltonian is in 
remarkable agreement with a full open-system calculation, even though the free-space geometry is 
far from the strong coupling regime of cavity QED. We explain this result using a simple model in 
which the atom couples weakly to the laser while coupling strongly to the vacuum. Additionally we 
place an upper bound on the total entanglement between the atom and all paraxial modes using a 
quantum trajectories unravelling. This upper bound provides a benchmark for atom-laser coupling. 

PACS numbers: 03.67.Mn,42.50.Ct,03.65.Ta,03.65.Yz 



I. INTRODUCTION 

The electromagnetic field describing a laser pulse is in- 
trinsically quantum mechanical and may become entan- 
gled with an atom when the two systems interact. The 
consequences of this fact for experiments which employ 
lasers as control fields has been the topic of several re- 
cent papers 0,0,0. All agree that these quantum effects 
are negligible compared to errors induced by spontaneous 
emission out of the beam. In particular, in our previous 
work 0, we calculated the rate at which information is 
gained about a two-level atom when one measures a near 
resonance laser pulse that has interacted with it. This 
measurement strength could be interpreted as the rate at 
which photons are spontaneously emitted into the modes 
defined by the laser beam - a small fraction of the total 
emission rate. Indeed, as commented by Itano Q, the 
quantum fluctuations in a laser beam are equivalent to 
those in the vacuum, and thus using the standard unitary 
Mollow transformation [5| , all quantum effects related to 
the atom-laser interaction in free space are derived from 
spontaneous emission in the transformed frame. 

Though not an important source of error in quan- 
tum operations, entanglement is nonetheless generated 
between the atom and the laser beam. Photons sponta- 
neously emitted into the laser modes change the quan- 
tum statistics of those modes in a way that is correlated 
with the state of the atom. This entanglement might 
prove useful as a resource in some quantum communica- 
tion or quantum computation scheme jfj. Moreover, by 
determing how much entanglement one can generate in 
free space, we obtain a benchmark for comparison with 
other entanglement generation schemes. For example, 
the entanglement generated in a cavity QED geometry 
could be compared to the entanglement that would 
be generated in the absence of the resonator. 

Quantifying the entanglement generated between an 
atom and a laser in free space is complicated by the mul- 
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timode nature of the field in a traveling wave geometry. 
The Hilbert space describing the laser modes is the tensor 
product of many independent subsystems. Additionally, 
the laser and atom do not form a closed system, as the 
other vacuum modes of the field interact with the atom. 
Since spontaneous emission into these unoccupied modes 
dominates over quantum features within the laser modes, 
it is generally inconsistent to neglect their effects when 
calculating atom-laser entanglement. 

To analyze this problem, we use the formalism devel- 
oped in our previous paper 0. We review the salient 
features in the next section. We then proceed to exam- 
ine the entanglement generated between the laser pulse 
mode and the atom. This result is compared to simi- 
lar results obtained using a closed-system single mode 
picture. We follow this with a calculation of an upper 
bound on the total entanglement generated between the 
atom and all modes propagating with the laser beam. A 
brief summary of results is presented in the concluding 
section. 



II. PARAXIAL FORMALISM 

In many problems it is useful to have a local description 
of the quantized field. For example, when the field inter- 
acts with a well localized atom or when an experiment 
uses a detector to continuously measure the field, a de- 
localized plane wave description is awkward. In systems 
where a physical quantization volume is defined, such as 
in cavity QED Q] , a local set of quantized modes is nat- 
urally defined. However, in free space no such intrinsic 
volume is present. A quantization procedure which con- 
siders local temporal modes rather than delocalized spa- 
tial modes has a limited regime of validity. Instead, 
through an appropriate coarse-graining of the standard 
modal decomposition, we can achieve the desired local 
picture |3(. A brief review of the pertinent points fol- 
lows. 

The procedure begins by separating the field modes 
into paraxial and nonparaxial modes. In the canonical 
example a laser beam impinges upon an atom, interacts 
with it, then illuminates a detector. This optical sys- 
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tem defines a set of modes traveling essentially in one 
direction with a small diffraction angle. We quantize 
these paraxial modes by associating a Fock space to every 
member of a coarse-grained set of wave packets of width 
Az and transverse area A. These form an approximate 
complete set for the paraxial subspace if the dynamics 
on a time scale smaller than At = c/Az are negligible. 
The coarse-graining length is arbitrarily chosen for con- 
venience and generally should not appear in any exper- 
imentally relevant result. It may be chosen to be the 
width of a trapped atom or to correspond to the detector 
bandwidth Aw, using Az = c/Aw. This latter definition 
is natural when considering the entanglement resource 
than can be detected. We introduce coarse-grained spa- 
tially localized mode functions 4>i{z) and an associated 
set of annihilation operators d;. The Hamiltonian for the 
paraxial field and atom in the interaction picture and 
under the rotating wave approximation is, 



H A F(t) = hg^20i(t) +a\a 



(1) 



where hg = dE vac = d 



AA2 



Here Oi(t) are window functions selecting out the mode 
currently interacting with the atom 



9i(f) 



(i - l)At < t < iAt 
otherwise 



(2) 



In this picture, each wave packet mode propagates 
past the atom, interacts with it according to a Jaynes- 
Cummings Hamiltonian [jj, then flies away. These dy- 
namics are intrinsically irreversible, in contrast to a cav- 
ity QED scenario in which excitaitons oscillate between 
the field and atom reversibly [?]. The information con- 
tained in these modes can act to decohere the atom at 
a rate re ~ Ta c g/A, where T is the total spontaneous 
emission rate and <j c s = 3tt/2}cq is the effective scat- 
tering cross section on resonance. The "measurement 
strength", re, represents the rate at which photons are 
spontaneously emitted into the paraxial modes. 

The nonparaxial modes which constitute the rest of the 
field are quantized using the standard free space quanti- 
zation procedure, associating a Fock space to every plane 
wave. In the following, these modes will be treated as 
a zero temperature reservoir into which the atom can 
decay. The only difference from standard spontaneous 
emission will be that the rate of decay is slightly reduced 
by the removal of the paraxial modes. 
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with Q p the solid angle subtended by the paraxial beam. 
Including the nonparaxial modes, the full master equa- 
tion for the atom-paraxial field system becomes 



dp 

~dt 



H A f,P 



V 



{o-+o--,/5} + r'<7_/5o- + , (4) 



where Haf is the Hamiltonian that couples the paraxial 
field modes to the atom, given in Eq. Q. 

In our canonical example the paraxial modes are taken 
to be in an initial coherent state, corresponding to the 
quantum state of a laser beam. We perform a Mollow 
transformation on the above master equation to re- 
move the initial coherent field amplitude of the laser and 
transfer it into the Hamiltonian. In the "Mollow pic- 
ture" , the initial state of the field is the vacuum and the 
master equation becomes 



dp —i 

~dt~ ~h 



1 r ; 

Hcdh + HxF,P — {<7+<X_,/5} + r <T_/5(T-| 



(5) 

where H co h is the coherent semi-classical Rabi-flopping 
Hamiltonian 



H, 



coh 



hga ((7_| 



(6) 



The Mollow transformation is achieved by a local unitary 
operator acting on the field alone and so will not affect 
the atom-field entanglement. Equation (J5J thus serves 
as our of basic dynamical description for entanglement 
calculations. 



III. QUANTIFYING ENTANGLEMENT 

In reviewing the paraxial formalism in the last sec- 
tion we made a distinction between the paraxial and 
nonparaxial field modes, the paraxial modes being those 
modes which travel past the atom and impinge upon the 
detector. We then traced out the nonparaxial modes to 
recover the master equation (JSJ) . The remaining paraxial 
subsystem is still a very large dimensional Hilbert space, 
composed of a tensor product of Fock spaces for each 
coarse-grained wavepacket in the pulse. Such a large 
Hilbert space has many subsystem decompositions, any 
one of which could be of interest when calculating entan- 
glement (Fig.^J. In this section we consider two possible 
decompositions: (i) entanglement between the atom and 
the mode defined by the laser pulse, (ii) entanglement be- 
tween the atom and the entire paraxial subsystem. The 
former is natural to consider when treating the laser pulse 
as a single mode, analogous to a Jaynes-Cummings cav- 
ity QED geometery pj ■ The latter is of interest when re- 
lating the atom-laser entanglement to the measurement 
strength of the laser and to the error in quantum logic 
induced by a control pulse. In addition, when calculated 
using this second division, we can find an upper bound on 
the atom-laser entanglement that can be used as resource 
for quantum information processing. 

To quantify entanglement, we use a monotone known 
as the "tangle" which can be related to the square of the 
more familiar concurrence and also to the entanglement 
of formation 0, ^J. For overall pure states where at 
least one of the systems has only two levels (is a qubit) 
the tangle is equal to the normalized linear entropy (or 
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FIG. 1: Two different modal decompositions of the paraxial 
field. The top line represents the local coarse-grained decom- 
position. The second two lines depict the Fourier decomposi- 
tion, which consists of many delocalized longitudinal modes. 
Only the symmetric mode (laser mode) is initially excited, 
starting in a coherent state, and is shown here shaded. The 
other nonsymmetric modes (two of which are depicted) are 
initially in the vacuum state. 

purity) of the marginal state of either subsystem |T^ , 

r(|^)(^|) = 2(l-T^[Tr B [|^)(^|] 2 ]). (7) 

This formula can be generalized to mixed states by taking 
the convex roof ^| 

T (p) = , m } n y^PiTQ^i) (V'il) P = y^Pi \i>i) O0i| • 

(8) 

We choose to use the tangle instead of the full Von Neum- 
man entropy due to its ease of calculation. 



IV. SYMMETRIC MODE ENTANGLEMENT 
A. Calculation 

We consider first the tangle between the atom an the 
"symmetric" paraxial mode (Fig. [I}, defined by the an- 
nihilation operator, 



a+ 



1 N 

— a 



(9) 



This mode defines a laser pulse of duration r. In calculat- 
ing the tangle in this decomposition all of the other field 
modes, nonparaxial and nonsymmetric paraxial modes, 
act as a reservoir which, when traced over, lead to a 
mixed-state description of our bipartite system. We as- 
sume kt <C 1, implying that the probability for the 
atom to emit a photon into one of the paraxial modes 
is very small. In the Mollow picture (used throughout) 
the Hilbert space of the quantum field may thus be trun- 
cated, allowing at most one excitation in the paraxial 
subsystem. Note that the atom is still allowed to spon- 
taneously emit an arbitrary number of photons into the 
nonparaxial field modes. Under these assumptions the 
atom and the symmetric field mode both behave as two- 
level systems, coupled together in an overall mixed state. 
The tangle for such two-level systems can be calculated 
using Wootters' general formula for two qubits 01 ■ 



We begin by calculating the atom-symmetric mode 
density matrix from Eq. (J5J for kt< 1. In principle, one 
could solve for the state of the entire atomic-paraxial sys- 
tem, and then trace out the nonsymmetric modes. How- 
ever, even limiting all paraxial modes to at most one 
total excitation there are still N + 1 possible field states, 
N being the number of coarse-grained modes. A more 
sensible approach is to keep only the "dynamic symmet- 
ric mode" at any instant (as defined below), tracing out 
any non-symmetric modes along the way. 

Consider first a single time step At = Az/c, with Az 
the coarse-graining length. At the start of this interval 
the system Hilbert space has dimension 2 ® 2. The sys- 
tem consists of the two-level atom and a quantized sym- 
metric mode of the field composed only of those coarse- 
grained modes which have already passed the atom, and 
restricted to one possible excitation. Assuming n < N 
modes have passed, the annihilation operator for this 
symmetric field mode will be 



i=i 



(10) 



After evolution over a time interval At, we choose a new 
subsystem division so that we again have a 2 £g) 2 system. 
The second subsystem will now correspond to the sym- 
metric mode over the n + 1 coarse-grained modes that 
have now passed the atom. The procedure may then be 
repeated to build up the full evolution, and the tangle 
can be calculated at each step. 

The procedure for evolving the state from time t to 
t + At is as follows. Suppose that n < N time intervals 
have passed. The density operator p n will be spanned by 
the basis states 



|0) + ,„ |5)|0) 
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I<7)|1) 



(11) 



During the next interval a new coarse-grained mode im- 
pinges upon the atom. This mode starts in the vacuum 
(in the Mollow picture) so the new state for this compos- 
ite system is 



P = Pn® |0)„ +1 (0| 



(12) 



Through the next time interval this state evolves under a 
combination of the Jaynes-Cummings Hamiltonian which 
couples the atom to the quantized mode, a semi-classical 
coherent atom-laser interaction, and spontaneous emis- 
sion into the nonparaxial modes. We model this by a 
short time integration of the master equation, Eq. JSJ. 
During this time the density matrix must be expanded 
to include the newly interacting mode. The basis states 
include the previous set, Eq. Qll|). tensored with |0) n+1 
for the newly added mode, plus the two new states 



10} 



u: 



n+l 



1.9) |0>_ 



1 



n+l - 



(13) 



States of the form |l) +n |l) n +i ar e higher order terms 
that are ignored under the assumption that kt< 1. This 
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gives a total of six basis states for the density operator. 
After performing the evolution over this time interval 
we trace out the nonsymmetric modes of the extended 
paraxial subsytem. To do this, we perform a unitary 
transformation on the field operators to create a new ten- 
sor product decomposition. Instead of the decomposition 
described in Eqs. i|llll3|) consisting of the n th symmetric 
mode £g> the (n + l) th coarse-grained mode, we choose 
symmetric and antisymmetric modes described by the 
annihilation operators 



1 



1 



(14) 



The transformation which achieves this is 



U = 




\° -a/5 



(15) 



Finally p n +i can be obtained by tracing out the antisym- 
metric field mode. 

The mixed state entanglement between the atom and 
the symmetric field mode can be calculated using Woot- 
ters formula for the tangle [10( directly from this den- 
sity matrix. The results will depend on several physical 
parameters: the area of the beam, the intensity of the 
laser \a\ 2 , and the duration of the pulse. Additionally, 
the numerical accuracy will depend upon the length of 
the time step used in the calculations, or equivalently 
the number of time steps N . To maintain the paraxial 
approximation, we choose the area of the beam to be 
A = A/a c f[ = 1000. This corresponds roughly to the 
quadrupole transition S-yo to D5/2 in 40 Ca + considered 
by van Enk and Kimble [l| , for which a laser focused to a 
spot size of 100 /Ltm 2 corresponds to A w 1600. For sim- 
plicity we consider only 7r/2 pulses, such that gar = tt/2. 
With the relations k — T/A and g = y/n/At — y/ kN/t, 
assuming a fixed N, we have a direct relationship between 
a and r. We choose as our free parameter the emission 
probability, IY, which sets the degree of entanglement 
between the atom and laser. 

A plot of the tangle versus Tt is given in Fig. [21 Each 
point in the plot corresponds to the total tangle measured 
just after the n/2 pulse has completely passed the atom. 
For very short pulses, the growth in the tangle is linear 
with Tt for all initial conditions, leading to a unit slope 
on the log-log plot. However the constant of proportion- 
ality is different. This difference can be attributed to the 
probability of the atom being in the excited state during 
the pulse interval - the more excitation the higher the 
rate of growth. With a greater probability of spontaneous 
emission, the behavior deviates strongly from this linear 
trend, peaking roughly where the pulse length and deco- 
herence time are comparable. After this time the entan- 
glement drops precipitously, and for at least one initial 
condition, the atom and symmetric field mode become 
separable. The maximum emission parameter shown in 




FIG. 2: The entanglement (tangle) generated between a two- 
level atom and a quantized tt/2 laser-pulse versus the deco- 
herence (IY) on a log-log plot. The pulse consists of N = 10 4 
coarse-grained modes and a normalized beam area A = 10 3 . 
Including spontaneous emission out of the beam, the tangle 
grows as a simple power law for IV <C 1, then drops expo- 
nentially fast for large decoherence. Shown for comparison 
are the results of the tangle calculated using a closed system 
Jaynes-Cummings dynamic that grown linearly for all times 
on this plot, but agree very well with the full calculation when 
the probability of spontaneous emission is small. Six different 
initial atomic conditions are shown: |e) (blue, solid cross), \g) 
(green, dashed), |e) +i \g) (cyan, solid), |e) —i \g) (black, dot- 
ted), |e) ± \g) (red, dash-dot). These last two initial condition 
show exactly the same tangle for all times (see text). 



the plot is Tt = vA = V1000, which corresponds to 
kt = 1/yA = 1/VlOOO; for a larger emission proba- 
bility, our truncation of the field Hilbert space to one 
excitation breaks down. 



B. Approximations 



The calculation in Sec. IIV Al yields the entanglement 
between the atom and the laser mode as it propagates 
past the atom. It is a complex procedure because we 
have explicitly separated the two sources of decoherence: 
spontaneous emission into the nonparaxial modes and 
spontaneous emission into the paraxial but nonsymmet- 
ric modes. We did this because spontaneous emission into 
the symmetric-paraxial mode is the source of the entan- 
glement and we must be careful not neglect this impor- 
tant effect. Nonetheless, emission into the nonsymmetric- 
paraxial modes is a small fraction of the total loss since 
s « f. We can consistently maintain the coherent cou- 
pling between the atom and the symmetric mode while 
neglecting emission into the nonsymmetric mode since it 
acts only as a small source a decoherence, a process dom- 
inated by off-axis emission. We thus lump all sources of 
decoherence together, with total rate T = k + T' . In this 
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case the Master equation becomes, 

dp % r 

— = — [H',p] + ~(2a_pa + -a + a_p- pa^a + ), (16) 

where H' = — h£l(a + + <5"-)/2 — hg(a& + + a)a_) is the 
driven Jaynes-Cummings Hamiltonian. This equation 
describes a cavity QED dynamic plus decay into a reser- 
voir When the evolution of the tangle is calculated 
using this equation of motion, the results are in complete 
agreement with those obtained in Sec. IIV Al to within 
numerical accuracy. 

A more radical simplification is to ignore the decoher- 
ence term in Eq. (|16fl altogether when the interaction 
time is small compared to the atomic lifetime. The re- 
maining dynamics represents a closed system, with the 
atom and the symmetric mode evolving unitarily. 

For this evolution, the initial bipartite pure state re- 
mains pure and the entanglement can be calculated in a 
straightforward manner. Such an approximation is rou- 
tinely made for atoms in cavities in the strong coupling 
regime, where g 3> T |7fl. In free space we can express 
g in terms of k; the strong coupling approximation re- 
quires k S> r 2 r, or equivalently TtA <C 1. The beam 
area may be chosen arbitrarily large, however, so that 
this inequality need not hold. One is not in the strong 
coupling regime in free space and so agreement between 
the entanglement calculated using closed and open sys- 
tem models should not be expected apriori. This has 
been the source of some recent controversy 0, 0, ^| . 

To compare the numerical results obtained using Eq. 
(J5J with those obtained using the closed system evolution 
Eq. I|17|) we first calculate the tangle analytically using 
the closed system dynamics. Though one could measure 
the closed system entanglement using the von Neumann 
entropy of the marginal state as was done in pj , for direct 
comparison with our results, we use the tangle given in 
Eq. (7J. A simple calculation yields, 

r = /tT (2K) +7 rf + (^(4-^) +Q 

where the Pauli matricies describe the two-level atom in 
the usual way. Expectation values in this equation are 
to be taken using the initial state of the atom. The re- 
sults of this approximation are shown in Fig. [21 For the 
parameters where all approximations are valid, we see a 
simple linear growth of the entanglement with emission 
probability. Though slightly larger, the entanglement 
calculated via the closed system dynamics is in very good 
agreement, 0(Tt), with our numerical results for the full 
atom-symmetric field mode tangle when Ft <C 1, as seen 
in Fig. This agreement holds even when the system is 
not in the strong coupling regime, ATt > 1. 

To understand why the entanglement calculated using 
closed and open system dynamics agree so well even when 
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FIG. 3: The decoherence channels for the atom-laser mode 
interaction are shown, with each block representing a subsys- 
tem of the entire Hilbert space. Subsystem A (the atom) is 
connected directly to reservoir R (all vacuum modes) with 
coupling rate T (the spontaneous emission rate). System B 
(the laser mode) is connected to the reservoir 7? only indi- 
rectly through system A. Since B couples to A at rate k (the 
measurement strength), it is indirectly coupled to R at rate 
Tut, where r is the duration of the interaction. 



the strong coupling assumption is violated, consider the 
processes which lead to the loss of entanglement. The 
inherent dynamics of the atom-laser mode system can be 
represented by a bipartite system, with subsystems A, B, 
coupled to a reservoir of other electromagnetic modes, 
R, as depicted in Fig. [31 System A is coupled to system 
B with strength k, while it is coupled to the reservoir 
R with strength T. System B is not coupled directly 
to the reservoir, coupling only indirectly through system 
A. The total system consisting of both subsystems A, B 
will start in some pure state. From the perspective of 
quantum trajectories 0], there are two possible ways 
in which the reservoir R can reduce the entanglement 
generated by the dynamics between A and B. In the 
first method the nonunitary but deterministic no-jump 
evolution of system A leads to a reduction of this en- 
tanglement with probability 1 — e~ rT . The total loss of 
entanglement is thus proportional to TtT for Tt <C 1, 
where T is the tangle generated in the closed system dy- 
namics between A and B. Since entanglement will be 
generated in the closed system at a rate proportional to 
the coupling k, the total loss of entanglement is of order 
0(Ttkt). Alternatively, a nondeterministic jump in a 
quantum trajectory can reduce entanglement if it leads 
to a particular statistical mixture of states as described 
below. The probability of this particular type of corre- 
lated quantum jump is, however, limited by the smaller 
decoherence rate of the two subsystems A, B. In this 
case system B can only decay through system A in the 
interval r at a rate proportional to Ttk. The amount of 
entanglement that can be destroyed through this process 
will then be of order 0(Ttkt). Thus both deterministic 
decay and quantum jumps lead to entanglement loss of 
order Ttkt. This is sufficiently small to ensure that the 
closed system dynamics generates the same entanglement 
as the open system dynamics when Tt <C 1 and kt < 1, 
independent of the relation between them, i.e. even when 
the strong coupling approximation breaks down. 
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The above explanation can be made rigorous for the 
atom-symmetric field mode system. Assume that Ft <C 1 
and kt <C 1, though ATt may be large. For Tt <C 1 the 
state of the system after the pulse has passed may be 
represented by a density operator of the form 

p = (1 - e) (V>| + ep jump . (19) 

Here \ip) is the state evolved under the closed system 
unitary dynamics, while /0j ump represents the state con- 
ditioned on a jump having occurred during the inter- 
action time. The probability for a jump is denoted 
e and will be proportional to Tt. These assumptions 
alone are insufficient to ensure that entanglement is pre- 
served. As an example, consider the entangled state 
|V>) = (1 - £/2) \gl) + |e0). This state has tangle T w 
4£ for £ -C 1. Note that the state is equivalent under 
local unitaries to any other state with the same entangle- 
ment. If we choose /3j urop — 1/2 |el) (el | + 1/2 \g0) (g0\, 
the statistical mixture p given in Eq. ljlH[) is separable for 
VT < e, as can be seen by its positive partial transpose 
|l7l Il8|. For our particular case, an atom in free space, 
V£ < 5 t j since entanglement is generated by the atom- 
field coupling. So for this choice of /5j ump , the system can 
be made separable when g < T, or equivalently ATt > 1. 

Since our calculations shows that the open system 
maintains the same entanglement as the closed system, 
even when ATt > 1, the dynamics must be restricted 
such that these special statistical mixtures are rare. To 
see this, note that the post-jump state we considered 
above must be generated by a highly correlated jump 
process. The probability of such a process is very low. 
Indeed, to good approximation, we expect /0j ump to be of 
the form 

Pjump = Patom <8> |vac) (vac| (20) 

where /5 at0 m is some arbitrary (possibly mixed) state of 
the atomic subsystem. This ignores the process of spon- 
taneous emission into a nonparaxial mode followed by 
subsequent excitation of a paraxial mode. This higher 
order process may be safely ignored when we are consid- 
ering entanglement of order kt, and ^/kt 3> ktTt, which 
is satisfied when Tt <C 1. Inserting Eq. lf2*U|) into Eq. ijl§|) 
the tangle may then be calculated using Wootters' for- 
mula, assuming that \ip) is an arbitrary state that can 
be generated using the closed system dynamic. To low- 
est order, O(ktTt), this agrees with the tangle obtained 
using closed system dynamics. 

C. Bounding the Total Entanglement 

We now consider the entanglement that is generated 
between the atom and all paraxial modes. This is the 
physically relevant quantity that determines the error 
rate in atom-laser control, modulo spontaneous emission 
into the nonparaxial modes. In addition, it gives a mea- 
sure of the dynamically generated entanglement resource 



that is accessible in principle. A calculation of this to- 
tal entanglement is significantly more complex than the 
calculation performed in Sec. II VI as the Hilbert space as- 
sociated with all coarse-grained paraxial modes is much 
larger than that of the single symmetric mode. Even if 
we assume kt « 1, and thereby restrict paraxial field 
modes to a single excitation as in Sec. IIVI the remaining 
Hilbert space has dimension 2 ® (N +1). A calculable 
measure of entanglement for a mixed state of a 2 ® N 
system for N > 2 is not known except for some special 
cases which do not apply here 0, |2(j . We thus resort 
to calculating a bound for the total entanglement, using 
the tangle as our measure. 

The general difficulty with calculating the entangle- 
ment of a mixed state arises because the ensemble decom- 
position of a density operator into a statistical mixtures 
of pure states p = J2iPi \4>i) (V'il i s n °t unique. Averages 
of a pure state entanglement monotone for different en- 
semble decompositions do not usually lead to the same 
result and one must thus perform a difficult minimization 
search to find the actual entanglement. Since the tangle 
is a convex function, however, any ensemble decomposi- 
tion provides an upper bound according to, 

T( P ) < 5>T(|^> (HM^Pi} (21) 

i 

Given a known ensemble decomposition, we can calculate 
the right hand side by averaging the tangle for the given 
bipartite pure states, over the appropriate probability 
distribution, using Eq. J7J). 

A natural ensemble decomposition is obtained by inte- 
grating the master equation of the atom-par axial system 
using the quantum trajectory method [16j. A stochastic 
Schrodinger equation describes the evolution of an open 
quantum system with nonunitary evolution due to its 
coupling to the environment, and conditioned on a ficti- 
tious measurement performed on the environment with 
unit efficiency. For each measurement record M(t) there 
is a corresponding pure state trajectory \ipM(t)}- Aver- 
aging over the all possible measurement records results 
in a mixed state description, with density operator given 
by the convex combination 

P(t) = EpmW |V>m(*)> (Vm(*)| • (22) 

M 

From this construction one may immediately write down 
an upper bound on the tangle using Eq. (|2*T)l . 

Brute force calculation of all the required trajectories 
is straight forward but prohibitive. If we choose the ficti- 
tious measurement to be photodetection, a measurement 
record will consist a sequence of bits denoting count/no- 
count decisions in each interval (t, t + At) (ignoring the 
possibility of two or more photodetections as it is of or- 
der St 2 ). Assuming there are N = 1000 such intervals 
then there are 2 1000 trajectories. Our calculation is made 
tractable, however, by noting that only the most recent 
jump in any trajectory will contribute to the tangle; the 
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state of any field modes which passed the atom prior to 
the last jump will not be entangled with the atom since 
the post-jump state is separable (the atom in the ground 
state and the field in a coherent state). The evolution af- 
ter a jump has occurred will be independent of when the 
jump occurred as the system is Markov, and all of the 
quantized paraxial modes prior to interaction with the 
atom are in the same coherent state. It is thus sufficient 
to calculate two no-jump trajectories - the initial trajec- 
tory \tp A (i)} an d the post-jump trajectory \ip B {i)). Both 
represent deterministic evolution according to a nonHcr- 
mitian Hamiltonian, including the coupling of the system 
to the nonparaxial modes assuming no quantum jumps 
into the nonparaxial modes. They differ only in their 
initial conditions: \rp A (i)) is the evolution given the ini- 
tial state of the atom at t = and is the evolution 
starting with the atom in the ground state right after a 
jump. An upper bound on the tangle then follows from 
Eq. (EU), 

N 

T{p) < Plj(i)T P j(N + P NJ T NJ , (23) 
1=1 

where -Plj(*) and Pjsij are respectively the probabilitiy 
that the last jump occurred during the time interval i and 
the probability that no jump occurred during the entire 
evolution. The tangle associated with no jumps during 
any interval follows from Eq. Q, T NJ = T(\^ A (N))), 
whereas the post-jump tangle is given by Tp j{N — i) = 

T(\MN -*)))• 

The relevant probabilities that appear in Eq. I)23|) may 

also be calculated using the fiducial trajectories. In par- 
ticular under nonunitary evolution without normaliza- 
tion, the probability that no jump occurs in a given tra- 
jectory up to time t is HV^t)!! 2 - Thus, the probability 
that no jump occurs during any of the N intervals is 
Pnj = IIV'a(-^)|| 2 - The probability that the last jump 
occurred in interval i is given by, 

^j(i) = P N AN\i)Pj(i), (24) 

where Pnj(N\i) = \\iJ)b(N— i)\\ 2 is the conditional prob- 
ability that no jumps occur between the interval i and 
N given that a jump occurred in interval i. The total 
probability that a jump occurred at time step i, Pj(i), 
will satisfy the following equations: 

Pj{i) = P A (i) + £ Pj{i\j)P A (j), (25a) 

3 

Pj(i\j) = Pb{* PMk)P B {k - j). (25b) 

k>j 

The first equation expresses the total probability for a 
jump in interval i as the sum of the probability that the 
system had its first jump in interval i, P A (i), plus the 
probability that it first jumped in interval j < i then 
some time later jumped in interval i. Pj(i\j) is thus the 
conditional probability for a jump at i given a previous 
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FIG. 4: An upper bound on the the entanglement (tangle) 
generated between the atom and all coarse-grained paraxial 
field modes is plotted versus the total decoherence (Ft) in a 
log- log plot. The different curves represent different atomic 
intial conditions with parameters given in Fig. |2] 

jump occurred in interval j. Equation (|25bfl provides 
a recursive relation for Pj(i\j) in terms of, Ps(i — j), 
the probability for a jump to occur i — j intervals after 
a previous jump occurred, with no intervening jumps. 
The probability P A (i) can be easily calculated from the 
fiducial trajectory \ip A (i)) and similarly Pp(i — j) can be 
calculated from \ips(i — j)) 01- Since Pj(i\j) is only a 
function of i — j, as follows from the Markov property, 
the recursion relation l|25b|l can be solved independently. 
We find a suitable solution by truncating the recursion 
at n steps corresponding to at most n + 1 jumps into 
the nonparaxial modes. Solving for PLjii), using this 
truncated solution, we can then obtain an upper bound 
on the tangle according to Eq. 

Our upper bound on the tangle is plotted in Fig. ED 
as a function of time. For short times, linear behavior is 
exhibited, as in the previous calculations for the symmet- 
ric tangle. Note, however, that the rate of entanglement 
generation is always greater than the symmetric calcula- 
tion. For short pulse durations Tt < 10~ 2 there is only a 
5% difference in rates for the initial condition \g) +i |e), 
the state that is most excited during the interval. All of 
the plotted upper bounds turn over at roughly the same 
point, slightly past where the pulse duration equals the 
decoherence time, and then converge on a single limit for 
large decoherence probability. This behavior is an arti- 
fact our calculation - for long times all information about 
the initial state has been lost since all trajectories with 
appreciable probability have undergone multiple jumps. 

V. SUMMARY 

Entanglement is generated between a laser beam and 
an atom in free space due to the atom's spontaneous 
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emission of photons into the paraxial modes, superim- 
posed on the laser's coherent field amplitude. Through 
the use of judicius approximations, we were able to quan- 
tify this entanglement, a fundamentally hard problem for 
a large dimensional open quantum system. In particular 
we calculated the entanglement between the atom and 
the field mode defined by the laser pulse. Of particular 
interest is how our measure of a entanglement is repro- 
duced under further simplifying assumptions. We con- 
sidered two models: a lumped decoherence model and 
a closed-system model. In the former all sources of de- 
coherence, paraxial and nonparaxial, are lumped into a 
single decoherence term. This is in excellent agreement 
with the full calculation when the paraxial approxima- 
tion holds, I » 1. This is as expected, since the vast 
majority of the decoherence is due to the nonparaxial 
modes. More surprisingly, a closed system model which 
treats the quantized atom laser-pulse dynamics via a sin- 
gle mode Jaynes-Cummings Hamiltonian is in excellent 
agreement when Tt -C 1, even when the strong coupling 
assumption was violated. This result was unexpected 
since the naive picture would argue that, in the absence 
of the strong atom-mode coupling, ATt > 1, spontaneous 
emission into other modes should destroy entanglement. 
We have shown, however, that the particular form of the 
coupling between the vacuum modes and the atom-laser 
system in fact preserves the entanglement. Thus, there is 
no contradiction between the observation that all quan- 
tum effects in the atom-laser interaction are due to spon- 
taneous emission 141 and the fact that atom-laser entan- 



glement can be accurately modeled using a single mode 
Jaynes-Cummings dynamic Q, Q when Tt <C 1. Such 
agreement is ensured by the fact that spontaneous emis- 
sion is the only decoherence mechanism. 

Finally, we placed an upper bound on the total entan- 
glement that is generated between the atom and the full 
set of paraxial field modes. We used the quantum trajec- 
tory method to calculate a dynamical unravelling of the 
master equation, which provides a natural ensemble de- 
composition of the bipartite density matrix. Averaging 
the pure state tangle over this ensemble leads to an easily 
calculated upper bound. It remains to be seen how tight 
this bound is. Nonetheless, it provides a useful bench- 
mark on the total entanglement resource that is available 
for quantum information processing applications. 
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